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We discuss homogeneous and isotropic cosmological models driven by SU (2) gauge 
fields in the framework of Einstein gravity. There exists a Yang-Mills field configura- 
tion, parametrized by a single scalar function, which consists of parallel electric and 
magnetic fields and has the stress tensor mimicking an homogeneous and isotropic fluid. 
The unique SU (2) gauge theory with spontaneous symmetry breaking sharing the same 
property is the Yang-Mills coupled to the complex doublet Higgs, this exists only in the 
case of the closed universe. This model contains an intrinsic mechanism for inflation due 
to the Higgs potential. Our second goal is to show that a successful inflation can be 
achieved also within the pure Yang-Mills theory adding an appropriate theta-term. 
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1. Introduction 

A remarkable fact is that the SU{2) Yang-Mills (YM) field admits isotropic and 
homogeneous configurations parameterized by a single scalar function. This is some- 
what surprising since the energy-momentum tensor of a single vector field in flat 
spacetime always contains anisotropic stresses EiEj, BiBj so this feature is clearly 
impossible in the U{1) case without spatial averaging. In the SU{2) case one has 
three vector potentials A° producing electric and magnetic component fields Ef ^ Bf 
which make the stress-tensor diagonal and isotropic if Ef = ESf and Bf = BSf , 
indeed EiEj = E'^SfSj = E^Sij. The electric and magnetic vectors are parallel, so 
both field invariants, the scalar T — F°^F"'^'^/2 — 3{B^ — E^) and the pseudoscalar 

g ^ gM^-Ar^a^^a^/g ^ g^;^ ^^.^ nOU-ZCrO. 

Such an ansatz can be generalized to any of the FRW metrics 

ds^ =dt^ -a{tf [dx^ + Y.k{de^ + sin^ 0d(l)^ )], (1) 
closed (fc = 1, El = sinx), open (fc = —1, = sinhx) or spatially flat (k = 
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0, So = and in the matrix form reads ^ 

F = F'^T"- = f (T„ dtAdx + Tgllk dtAde + T^T.k sin 0dtA dcf) 

+I]fe(/2 - k) (T^ dxAde- Tg sin 9 dx A dcj) + r„Sfe sin 0d0A dcf) . (2) 

Here the rotating SU{2) generators are used: 

T„ = T'^n72i, Te = T'^e72z, - r''ey2i, (3) 

where n", e^, are the spherical unit vectors, and r" are Pauh matrices. This prop- 
erty remains vahd also for larger gauge groups containing an embedded SU(2) 
Remarkably, the EYM cosmology is also solvable at the quantum level, the quan- 
tum FRW cosmology having been discussed in a number of papers in the 90-ies 
An interesting feature of the EYM quantum cosmology is a possibility of tunnel- 
ing transitions between de Sitter and hot FRW cosmologies Similar models were 
formulated in the context of supergravities In the spirit of the string theory the 
Born-Infeld generalizations were also considered^. More recently, the YM cosmo- 
logical models with mo dified lagrangians were suggested in the context of the dark 
energy models I 9 | 10 | ll l. 

Now consider the gauge theory with spontaneous symmetry breaking. For the 
SU{2) group one has two basic possibilities for the coupled Yang-Mills-Higgs (YMH) 
system: the real triplet Higgs and the complex doublet Higgs. The non-interacting 
Higgs belonging to any representation of the gauge group can form isotropic and 
homogeneous configuration. But for an interacting YM-Higgs system it is not so. 
It turns out that no homogeneous and isotropic configuration of interacting YM- 
triple Higgs exist which would be compatible with the Friedman-Robertson- Walker 
metric for any value of k. I n the complex doublet case such a configuration exists 
in the closed case A; = 1 This has a simple explanation: The doublet Higgs has 
SO (4) symmetry which is the isometry group of the three-sphere. This is another 
interesting coincidence, which selects the YM-doublet Higgs as natural "realistic" 
field-theoretical model for cosmology. This model can accommodate for infiation, 
and at the same time naturally introduces a vector field which could play a role of 
curvaton 

2. Yang-Mills-Higgs cosmology 

Consider the EYMH action with complex doublet Higgs: 

S = j [--^R If;^F'^'- + liD,^)^D,^ ^ ($t$ - ^^4^^ 

where <i> is the complex doublet Higgs, and 1?^$ 9^$ + 5yl°Ta$. The YM matrix 
potential — A'^ljTa generating the field strength for fc = 1 can be chosen as 

AlTa^^—^Ud,U-\ C/ = e2^^", (5) 
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while an ansatz for Higgs is 

$ = /i(<)C/$o, H^o = 1, (6) 

where $0 is an arbitrary norniahzed constant spinor. The compatibihty is this con- 
figuration with the required spatial symmetries follows from the identity 

The EYMH action contains three different mass parameters: the Planck mass, the 
mass of the W-boson, and the Higgs mass 

Mpi = Mw^ gv, Mh - %/Ai;. (8) 

V G 



We rescale the Higgs scalar as h hMpi and introduce dimensionless parameters 



S'ym ^ 


- 5h + Sint, 




/^(- 




= 






J I 2N 



gMpi Mpi Mw g 

Finally, introducing the reduced Planck length I = l/{gMp\), we present the metric 
in dimensionless terms: 

ds^ {-N^dt^ + [dx^ + sinxHde^ + sin^ 0dip^)]} . (10) 

Substituting the ansatze (j6l5l) . integrating over angles and dividing by the volume 
of the unit the three-sphere, we obtain the one-dimensional reduced action S =^ 

S,r 

-—)dt, ^YM^^yi^^ — jdt, (11) 

!^{h^ - a^fNa^^ dt, 5int = -\ j '^^(Z + "^fNadt, (12) 

and the total derivative in the gravitational term is omitted. 
Variation with respect to N leads to the constraint equation 

a'^ J 6 a 
and we fix the gauge = 1 afterwards. For the energy density we obtain: 

e^\h'+^^ + Vj + Vh + Mnt, (14) 
where the potentials read: 

3(/2_l)2 p2 3;^2/y^;L)2 

The acceleration equation is obtained by variation of the action over a with account 
for the constraint: 

a 47r , 



^' + 72)=^^' H^-, (13) 
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leading to the following expression for pressure: 

P=lh' + i, + lVf-V,, (17) 



and therefore, 



a 6 \ Za 



h' + ^-Vk + Vf\. (18) 



Note that Vh and Vf enter with different signs, while the interaction term V^nt does 
not enter the acceleration equation at all. The field equations for the YM and Higgs 
scalar functions read 

3 
2^2' 



h + 3Hh = '^h(f + 1)2 - ^2/^2 _ ^2^^ 

3,9, „ 6 



f + Hf = -^h\f + l)-^{f-l)f. (19) 

It is worth noting, that two roots / = ±1 of the YM potential Vf correspond to 
two neighboring topological vacua. To show this, consider the Chern-Simons 3-form 



^3 = ^Tr [AAclA- —A A A A A] , (20) 



satisfying the equation 

dw3 = ^TrF A F. (21) 

Substituting here the potential ([5]) one finds that it is non-trivial in the closed 
universe k — 1, giving the winding number of the map SU{2) — > S^: 

Ncs^ f c.3 = i(/ + 1)^(2-/). (22) 

The vacuum / = — 1 is therefore topologically trivial: Ncs — 0, while the vacuum 
/ = 1 is the next non-trivial one with Nqs — 1- This does not mean, however, that 
both will be classical solutions to our system. 

One can see that the above equations have the following particular solutions: 

• If Higgs is in the false vacuum ft, = 0, we have the EYM system with the 
cosmological constant A = a^/32/4; the corresponding scale factor can be 
found in a closed formS21. The Euclidean sector of the same theory contains 
instanton and wormhole solutions interpolating between the de Sitter and 
the hot universes 13]. 

• A particular (unstable) analytic solution is cosmological sphaleron in the 
false Higgs vacuum: 



0, / = 0, a^^^. (23) 



This solution describes the YM field sitting at the top of the potential 
barrier separating two topologically neighboring YM vacua (note that f — 
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is not vacuum for YM ) and generalizes the pure EYM sphaleron of Ref . 
In this case H — and the universe is in the steady state. 

• There is another non-trivial steady state with non-zero /, h: 

• If YM field is in the trivial vacuum / = — 1, we have a pure Higgs k — I 
cosmology, which is the basis of the standard inflationary models. 

At the same time, there are several regimes which could be expected, but actually 
are not realized as classical solutions of the EYMH system: 

• The system never reduces to the pure EYM cosmology with the conformal 
equation of state p = and the hot universe metric as described in Ref. |T] 

• The non-trivial vacuum / = 1 is not a solution to our system: the interaction 
term does not vanish in this case, so the YM equation will not be satisfied. 
Thus, the YM field dynamically coupled to Higgs will have the unique 
vacuum state / = — 1. 

• The Higgs vacuum h = a ia not a solution unless the YM field is vacuum 
/ = -!• 

The system possesses the slow-roll regime when kinetic terms /, h are small relative 
to potential terms. There are three different slow-roll regimes: 

• infiation, with an approximate equation of state p = —e in the case of the 
dominant Higgs potential Vh; 

• hot universe with p = e/3 in the case of dominant YM "potential" V/; 

• zero-acceleration regime with the equation of state p — — e/3 in the case 
of dominant interaction term T^nt- The emergence of such an equation of 
state (usually referred as the string gas equation) regime is somewhat unex- 
pected, it also occurs in the high energy limit of the non-Abelian Born-Infeld 
cosmology 

Note, that Vf and T^nt are not really the potentials: they depend on the cosmological 
scale factor. 

The solutions ([23l) and ([M)) are fixed points of the six-dimensional dynamical 
system associated with the equations of motion. The first one corresponds to the 
local maximum of the Higgs potential, in which case interaction between the Higgs 
and the YM fields is switched off. The second point corresponds to the minimum 
of the potential /imin < a, shifted by interaction with YM. Correspondingly, the 
position of the YM function / is also shifted from the extrema of Vf . In this case 
the interaction plays crucial role. The analysis of perturbations around the fixed 
points shows that the first one is an unstable node. For the second fixed point one 
obtains stable and unstable directions in the phase space depending on the value of 
the parameter /3, see Ref. [12| for details. The most interesting case corresponds to 
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stability of this point against small perturbation, which translates in the existence 
of long quasistable states which terminates by an inflationary period. The number 
of e-folds during inflation depends on the initial value of the Higgs field in the 
quasistationary phase. 

Qualitatively, behavior of the integral curves obtained in numerical experiments 
is somewhat reminiscent of the hybrid inflation with two truly scalar fields. 
The interaction potential Vint plays the role of trapping potential which drives the 
inflaton to climb the Higgs potential. As a result, the YM field can substantially 
enhance the scalar inflation. 



3. Theta-inflation 

Consider now the pure YM theory described by the action S = J L^—gd^x, where 
the lagrangian L{T , Q) depends arbitrarily on two YM invariants 

T = -Fl^F''^''' 12, = -F"^''F;^/4, F"^'^ = gMjAr^^ ^^^^ 

We will focus here on the dependence of the lagrangian on the pseudoscalar Q. In 
gauge theories the linear in Q term is induced by instantons and is called theta-term. 
We will assume, however, more general dependence of L on motivated, e.g. by 
vacuum polarization. 

In the Refs. ITUl IT^ we derived the equation of state for FRW cosmologies driven 
by the YM field configuration ([2|) for a general lagrangian L(T , G)- The contribution 
coming from Q was found to produce the energy density and pressure with the ratio 
w = e/p=— Litis not difficult to show that this does not depend on the particular 
configuration of the YM field. The linear functional Sg = J Qy/—gd'^x does not 
depend on the metric: 

1 f e''"^^ ^ „ , 1 



SQ = -^j ^-^F^^FxrV^d'^x = ~ 2 y e'^'^^^F^.F^rd'^x (26) 

and hence does not contribute to the energy-momentum tensor. However, for more 
general dependence of the lagrangian on Q, it is not so, since 

^ = \Q9,u. (27) 
Using this, we obtain the YM energy-momentum tensor in the form 

where the second term looks like the variable cosmological constant. The field in- 
variants for the YM configuration ([2]) read 
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We also introduce the effective electric and magnetic fields 



J_ 



n = 



so that F = 3(£2 - H^), g = sen. 



(30) 



Then in the context of the FRW models we find for the energy density and pressure: 



'dQ 



P 



L 



dF dg 

Consider now the lagrangian with a "potential" depending on Q : 

T 
"2 

Then 3 , 



ViG). 



and the Friedman equations read: 



eg=V- 3£nV', 



(31) 



(32) 



(33) 



Stt 



3 

Stt 



{eg - Sjr) 



-=H + H- = 
a 6 

The dynamics of the YM field in this model is governed by the equation: 



d 
dt 



a 



df 



a 



df df 



(34) 



(35) 



Using the relation V" = —3^7, where primes denote derivatives with respect to Q, 
after some rearrangements we obtain: 



£ 



/ + 



Hf - 2/H 



1 + 



(36) 



We will investigate the regime when e'g terms dominates in the squared brackets, 
then the Eq. ([55]) is simplified to 

/-3i// = 2/^. (37) 

Contrary to the scalar inflation generated by the potential term, here the desired 
regime rather corresponds to the kinetic (electric) term being of the same order as 
the potential term (magnetic). Indeed, the theta-term, whose dominance produces 
inflation, is the product of the electric and magnetic components. The dynamics will 
depend on the initial state of the system, characterized by the ratio 'H/£, which 
should not be too small or too large. Let us investigate the slow-roll regime, which 
is expected if 



• Hubble friction term in the YM equation is dominant; 
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• the accelerating term in the Friedman equations is dominant: ejf <C eg; 

• the Hubble parameter is varying slowly, e = —H/H^ <C 1,7] = 
~H/{2HH) < 1; 

• the curvature term is negligible k/a^ ^ H^. 

In what follows we set fc = 0, just noting that the case fc = 1 is rather different and 
very interesting, it will be considered elsewhere. For the slow-roll approximation it 
may be not correct to use just / ^ 3iJ/. It will be better to work with the ratio 
ip = f /a, which is connected with the magnetic component via H — —ip"^ and whose 
derivative is related to the electric one: £=?/' + Hip. Thus we arrive at the system 



3 



ip + 2Hip = -tA 
(e^ + eg), H 



i/2 



3 



{eg - ejr) 



(38) 



The dominance of the theta-terms eg 3> ejr implies smallness of the first slow-roll 
parameter: 



< 1. 

ejr + eg 

To simplify further calculations we proceed with the potential 

vig) = -0ging, 



(39) 



(40) 



in which case the theta-energy density eg = 9Q is linear in Q, the same as the 
standard T term. Then the condition ej^ <^ eg of the accelerated expansion can be 
satisfied in a wide range of initial conditions and will persist long enough if 3> 1. 
In this case the condition e'g = 6* ^ 1 is satisfied and we arrive at the system (p8|) 
admitting an accelerated solution if the electric and magnetic components are of 
the same order. 

To ensure slow variation of the Hubble parameter H and the effective infiaton 
ip during the inflation one needs slowness of the parameter 5 - 
which the electric and the magnetic components read 



-j^, in terms of 



n 



= 25. 



(41) 



By virtue of the Friedman equations, H 



1 



V = 



e + S€ + 36^ 



(5) ' HH 

£jF, so we can rewrite the system (j38p as 

esn 



e = 2 



1 - 



H{l-6) 



2'H'-5 



(42) 



To ensure slow variation of (5, S/{HS) ^ 1, we set 6 
equation in the leading order in e reduces to 

e g2 _ 2^2 
2 £2 + 7^2 ' 



V 



-e/2 + 0(e2). Then the last 



(43) 
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therefore 77 ~ e, e ~ He^. Finally, smallness of e implies the constraint on the initial 
state and the parameter 9: 

^«'- (") 

so 9 should be large enough, with both electric and magnetic fileds non- vanishing. 
The number of e-folds is given by 

Using the relation between variations of ^ and H due to 5 = — e/2, we obtain 

V' [Hi , X 



where the initial values must satisfy the constraint 

- f - 8.^ff.^f , ^ (47) 



The slow roll ends when the condition pi)) fails. During the inflation the Hubble 
parameter decreases since H = —IGnej^/S, therefore the inflaton tp goes up. The 
electric component decreases as f ~ Hip ^ while the magnetic increases as 
\H\ - V^- The slow roll regime stops when the ratio H/S approaches the value 
9, i.e. the imbalance between the electric and the magnetic components becomes 
comparable with the imbalance of the standard and 0-term in the lagrangian. This 
happens when 

^--0 ^ ^l-8n9'iP!. (48) 

Then the dominance of the e'gE /H term at the right hand side of the Eq. p6l) ends, 
and the system dynamics changes. 

To estimate the number of e-folds, we substitute dH/H — ~2dip/i' into (|^5|) 
using the Eq. pSj) to relate 0e and (pi. For e use the system (|^ expressing it in 
terms of ip alone via During the main part of the inflation stage one has H > £, 
so that 1/e ~ -9£/{2'H) = 9H^ipf / {2ip^). Finally we obtain: 



Unlike the scalar case when the scalar field is rolling down, the YM inflaton ip 
is rolling up during the inflation. To avoid going beyond the Planck mass we 
should properly choose the initial conditions. Taking as the value of ipe the rescaled 
Planck mass gAfpi, we obtain from (|48p for the initial value in the same units: 
ipi = (87r6'^)^^/^ ^ 0e ^ 1- Hence the number of e-folds will be 

A^e ~ -(8716*2)3/5 2.3 61^/^ (50) 
3 

Therefore we need the parameter 9 to be just the order of ten to obtain the slow roll 
regime with sufficient number of e-folds. The strongest inflation takes place when 
electric and magnetic YM components initially are of the same order. 
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4. Outlook 

Incorporating gauge field into cosmological setting seems natural and promising. 
We have shown that homogeneity and isotropy selects the YM-doublet Higgs as 
the unique gauge theory with spontaneous symmetry breaking compatible with the 
closed FRW cosmology which admits both the inflationary and the hot universe 
regimes. Our second proposal consists in using a suitable theta-term to achieve 
cosmic acceleration with the pure YM lagrangian. 
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